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The Pomeranchuk instabihty, in which an isotropic Formi surface distorts and becomes anisotropic 
due to strong interactions, is a possible mechanism for the growing number of experimental systems 
which display transport properties that differ along the x and y axes. We show here that the gauge- 
gravity duality can be used to describe such an instability in fermionic systems. Our holographic 
model consists of fermions in a background which describes the causal propagation of a massive 
neutral spin-two field in an asymptotically AdS spacetime. The Fermi surfaces in the boundary the- 
ory distort spontaneously and become anisotropic once the neutral massive spin-two field develops 
a normalizable mode in the bulk. Analysis of the fermionic correlators reveals that the low-lying 
fermionic excitations are non-Fermi liquid-like both before and after the Fermi surface shape distor- 
tion. Further, the spectral weight along the Fermi surface is angularly dependent and can be made 
to vanish along certain directions. 



I. INTRODUCTION 



To a surprising extent, the electronic properties of 
most metals are well described by Landau Fermi liq- 
uid theory. The key building blocks of this theory are 
quasi-particles and a Fermi surface, the surface in mo- 
mentum space that demarcates filled from empty states. 
For a system of weakly interacting electrons, the shape 
of the Fermi surface is isotropic. On the other hand, 
if the interactions among the quasi-particles are strong 
enough, other possibilities might emerge. Pomeranchuk 
[1] showed that forward scattering among quasi-particles 
with non-trivial angular momentum and spin structure 
can lead to an instability of the ground state towards 
forming an anisotropic Fermi surface. For example, the 
simplest form of a two-dimensional nematic Fermi liquid 
can be thought of, in the continuum limit, as a realization 
of the Pomeranchuk instability in the angular momentum 
I = 2 and spin s = particle-hole channel whereby a cir- 
cular Fermi surface spontaneously distorts and becomes 
elliptical with a symmetry of a quadrupole (d-wave). 
Such an instability, depicted in Figure 1(a), is often re- 
ferred to in the literature as the quadrupolar Pomer- 
anchuk instability^. In practice, however, the instability 
breaks the discrete point group symmetry of the underly- 
ing lattice. As such, systems undergoing a Pomeranchuk 
instability exhibit manifestly distinct transport proper- 
ties along the different crystal axes. Experimentally sev- 
eral instances of anisotropic transport in strongly cor- 
related electronic systems have been reported [2-5] and 
occupy much of the current focus in correlated electron 
matter. Although it is unclear at present whether the 




^ The instabilities in other channels, ineluding the antisymmetric 
spin channel, are also interesting. (Note that there is no Pomer- 
anchuk instability in the I = 1, s = Q channel.) In this paper, 
however, we only focus on the quadrupolar Pomeranchuk insta- 
bility and the resulting nematic phase. 



FIG. 1. (a) Examples of the Pomeranchuk instability in (a) I = 2, 
(b) I = 3 and (c) i = 4 channels. The -|- and — signs on the 
anisotropic Fermi surfaces show the momentum gain and loss com- 
pared to the isotropic Fermi surfaces (which are shown above by 
the dashed lines). 



Pomeranchuk mechanism is the only root cause^, it nev- 
ertheless offers a framework where such anisotropics can 
be characterized in detail. So far, much of the work in 
this direction has focused on reaching the nematic phases 
of correlated electron matter, via the Pomeranchuk insta- 
bility, in a weakly interacting Fermi fluid [9, 10]. A nat- 
ural question to ask is whether a nematic phase could be 
approached, via a Pomeranchuk instability, from a non- 
Fermi liquid phase. To answer the question, one has to 
address the origin of the Pomeranchuk instability beyond 
the Fermi liquid framework. Unlike Fermi liquids, there is 
no notion of stable quasi-particles in non-Fermi liquids. 
Nevertheless, non-Fermi liquids have sharp Fermi sur- 
faces. Consequently, a Pomeranchuk instability, viewed 
simply as a spontaneous shape distortion of the underly- 
ing Fermi surface, is still a theoretical possibility. 

We explore here how such an instability can be stud- 
ied using the gauge-gravity duality, hereafter referred 
to as holography. The duality, which maps certain d- 
dimcnsional strongly coupled quantum field theories to 



•^ Nematic phases could also be obtained by the melting of stripe 
phases; see [7, 8] and references therein. 



{d + l)-dimcnsional scmiclassical gravitational theories 
in asymptotically AdS spacetimes, is a promising tool 
in modeling some features of strongly correlated con- 
densed matter systems. For reviews on the applications 
of holography for condensed matter physics, see [11-17]. 
We take the system to be d = 2 + 1 dimensional, as 
this is the most interesting case experimentally [2-5]. 
Also, wc consider the case where the broken phase is 
a nematic (non-)Fermi liquid with the symmetry of a 
quadrupole. In other words, we only explore the possi- 
bility of a quadrupolar Pomeranchuk instability. Such an 
instability is characterized by an order parameter which 
is a neutral symmetric traceless tensor Oij where i and j 
denote spatial indices. Thus, in order to use holography 
to realize the transition to such a phase in the boundary 
theory, one has to turn on a neutral massive spin-two 
field ifp^i, in the bulk and show that, under some circum- 
stances, it develops a non-trivial normalizable profile in 
the bulk radial direction. From the point of view of the 
physics in the bulk, there is a serious hurdle that one 
needs to overcome in order to successfully describe the 
propagation of a massive spin-two field, which has to do 
with the existence of extra unwanted degrees of freedom 
(ghosts) and superluminal modes. Eliminating ghosts 
and superluminal modes within a theory that also allows 
for the condensation of the operator dual to the massive 
spin-two field is the first step. In addition, a coupling of 
the correct type must be introduced in the bulk between 
the massive spin-two field and the fermions in order to 
bring about the shape distortion of the Fermi surface in 
the boundary. We present here a minimal model which 
can address each of these problem, and show that the 
isotropic Fermi surface of the boundary non-Fermi liquid 
gets spontaneously distorted and becomes elliptical. We 
also find that the broken (nematic) phase is a non-Fermi 
liquid. We interpret this transition as a holographic re- 
alization of the quadrupolar Pomeranchuk instability in 
non-Fcrmi liquids. Nonetheless, a disclaimer is neces- 
sary here. The background we study cannot be extended 
consistently to T = on the count that so doing would 
require the back reaction of the gauge field. Such a back- 
reacted metric would be non-Einstein and there is no 
causal Lagrangian for a spin-two field in a curved back- 
ground that is not Einstein. Within the background we 
consider, the vaccuum expectation value of the spin-two 
field actually grows as the temperature is increased. Con- 
sequently, there is strictly no instability in the traditional 
sense in the context of spontaneous symmetry breaking. 
Hence, our usage of the term instability is not entirely 
warranted. Nonetheless, the program wc outline here 
does provide a consistent holographic framework to ad- 
dress how to engineer shape distortions of the underlying 
Fermi surface. Hence, it serves as a valid first step in 
realizing the physics of the Pomeranchuk instability. 

The outline of this paper is as follows. Section II sets 
up the Lagrangian for the neutral massive spin-two field. 
Included here will be an explicit mechanism for obtaining 
a normalizable profile for the bulk spin-two field. Section 



III contains the fermonic degrees of freedom as well as the 
results illustrating that the low-lying excitations are non- 
Fermi liquid in nature as well as the shape distortion of 
the underlying Fermi surface. 



II. 



BACKGROUND 



Since we want to explore the possibility of a transition 
from a non-Fermi liquid to a nematic Fermi liquid phase 
in the boundary, the first step is to construct a back- 
ground which consistently incorporates a neutral mas- 
sive spin- two field (p^^. This field is dual to the operator 
whose vacuum expectation value (vev) is the order pa- 
rameter for the nematic phase. 



A. BGP Lagrangian 

In fiat {d+ l)-dimensional Minkowski spacetime, Fierz 
and Pauli suggested [18] a Lagrangian for a free neutral 
massive spin-two field which is quadratic in derivatives 
and correctly describes the propagating degrees of free- 
dom for such a field. Generalizing to curved spacetimes, 
Buchbinder, Gitman and Pershin [19], following [20], pro- 
posed a simple Lagrangian (quadratic in the massive 
spin-two field, and also quadratic in derivatives) which 
describes the causal propagation of the correct number 
of degrees of freedom of a neutral massive spin-two field 
in a fixed Einstein background in arbitrary dimensions. 
In four bulk dimensions, their Lagrangian for (Pfj,^, de- 
noted by >CbgP: reads 



^BGP 



+ 2V^(^,pV''(^'^^ - 2V^(^^,V V 
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Here R is the Ricci scalar. Also, we have defined ip — (f'f^. 
Several comments are warranted here. First, it is crucial 
that the fixed background satisfies the Einstein relation 
(which, in four bulk dimensions, reads) R^i, = Ag^i,, with 
A being the cosmological constant, otherwise some con- 
straint equations become dynamical which would give 
rise to extra propagating degrees of freedom for (^^^. 
In our context this restriction simply implies that one 
should ignore the backreaction of ipf^iy, as well as other 
bulk fields that will be introduced later, on the metric of 
the fixed Einstein background. Second, in the fiat space- 
time limit, the above Lagrangian becomes the Fierz-Pauli 
Lagrangian [18]. Third, the massless limit of the BGP 
Lagrangian, i.e. m? ~ 0, has an emergent gauge sym- 
metry (diffeomorphism) which is unique to the graviton: 
Lp^i, -^ <P/ii/+V(^^^) with ^^ being an infinitesimal vector. 
Indeed, as pointed out in [19], this symmetry has been 
used to settle the ambiguity associated with the definition 



of the mass of f^,^. Fourth, for m? = R/Q (here d+1 — 4) 
the above Lagrangian has another emergent gauge sym- 
metry [21]: ip^i, -^ ip^i, + {V^V^+gf^i,R/l2)e with e being 
an infinitesimal scalar gauge parameter. In the following, 
wc always assume that rn^ ^ and rr? ^ R/Q. Exclud- 
ing these two special values for the mass, the equations 
of motion and constraints obtained from (1) are given by 
[19] 
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The above expressions give the correct number of propa- 
gating degrees of freedom for a massive spin-two field in 
(3-|-l)-dimensions. Such a field transforms in the five- 
dimensional irreducible representation of S0(3), with 
S0(3) being the little group of the Lorentz group 
S0(3,l). 

As we alluded to earlier, in order to satisfy the con- 
straint R^i, = h.g^vi we will work in a regime where mat- 
ter fields do not backreact on the bulk geometry. So, 
we take the Schwarzschild AdS4 black hole as our back- 
ground geometry (A = —3) 



ds = g^i^dx^dx^ 

= r^ (-/(r) dt^ + dx^ + dy^) + 



dr^ 



r^f{r)' 



(7) 



where /(r) = 1 — (rp/r'') with tq being the horizon radius, 
given by the largest real root of /(ro) = 0. Note that we 
are working in units where we set the curvature radius L 
of AdS4 equal to unity. Also, note that in the coordinates 
we have chosen above, the asymptotic boundary of the 
spacetime is at r — >■ oo. The temperature of the black 
hole is given by T := 3ro/(47r). 

To solve for (pp,i,, we only consider the configuration 
where ip^i, — (pp,,v{r) for all n^v £ {t,r,x,y}. Analyzing 
the equations (2)-(6), one can show [22, 23] that it is 
consistent to put (^^^(r) in the form 



Vf^i^ir) 





tp,j{r) 



where ipij{r), with i,j — x,y, is a symmetric trace- 
less two-by-two matrix. Defining the "director" (f{r) — 
'^xxi'f) + "i-y^xyif), under a rotation 9 in the {x,y)-p\ane, 
ip{r) transforms as 



ip{r) -> e"*V(^)- 



(9) 



One can then choose a particular value of 9 to set either 
'^xxi'i') or fxyif) equal to zero, or make them equal. 



Suppose now there exists a (2-|-l)-dimensional field 
theory dual to the gravitational setup under consider- 
ation here. Focusing for the moment on the bulk neutral 
massive spin-two field (ff^,, holography tells us that it 
should be dual to some operator in the boundary the- 
ory. Since (p^i, is a bulk field with spin greater than zero, 
its number of components does not naively match the 
number of components of the dual operator. But, as ex- 
plained in [23], since ip^i, can be put in the form (8), it is 
(pij which sources a neutral symmetric traceless operator 
Oij in the boundary theory. (The vev of the operator 
Oij determines the nematicity of the (2-|-l)-dimensional 
boundary theory.) More concretely, equation (2) implies 
that (pij{r) takes the following form near the boundary 



^,j{r)^A,r^-^{l + ---)+B,^ 



„2-A 



(1 



(10) 



with Aij and Bij being symmetric traceless tensors. Aij 
is the source for the boundary theory operator Oij while 
Bij is proportional to its vev, (Oij). In (10), A denotes 
the UV scaling dimension of Oij , which is related to the 
mass of the bulk field (pfj_i, via 



2 



9 

4 



(11) 



Note that for m = 0, i.e. when ipf^i, is the graviton, equa- 
tion (11) yields A = 3 which is the dimension of T^^, 
the energy-momentum tensor operator, of the (2+1)- 
dimensional boundary theory. There is a Breitenlohner- 
Freedman bound [24] for the propagation of (p^i, in an 
asymptotically AdS4 geometry which reads rn^ > 
[23, 25, 26]. Equation (11) then implies that A > 3. 
In our discussions in this paper, we will always take the 
mass squared of ip^^i, in the asymptotic AdS4 region of 
the geometry to satisfy the condition m^ > 0. Note that 
the condition m^ ^ R/6 ~ —2 is then trivially satisfied. 

Given the ansatz (8), we are interested in solutions in 
which ipf^i,{r), or rather ipij{r), is regular near the horizon 
and normalizable as r ^- oo. More concretely, near the 
horizon, we look for a solution of the equation (2) which 
is regular as r — > rg , namely 



•Piji^) 



+ bi]{r - To) + c^j{r - r^f + 



(12) 



where the coefficient tensors bi 



are all deter- 



mined in terms of aij. In the boundary as r — > oo, we 



(8) demand fijir) to be normalizable 



^,,{r)^B,,r^-^ {! + ■■■) 



(13) 



If a solution with the above two boundary conditions ex- 
ists, then the boundary theory would be in a phase where 
the neutral symmetric traceless operator Oij (which is 
the operator dual to ipij) spontaneously condenses. It is 
easy to show that given the set of equations (2)-(6), such 
a solution does not exist. As we will explain in the next 
section, one way forward is to minimally modify the BGP 



Lagrangian (1) such that the new action yields the cor- 
rect number of causal propagating degrees of freedom for 
iPfiv, and also permits the aforementioned spontaneous 
condensation to occur. 



B. Modifying the BGP Lagrangian 

Undoubtedly, there are many ways to modify the BGP 
Lagrangian to facilitate spontaneous condensation of the 
operator dual to >Pfj,u- Instead of categorizing all such 
modifications, we take perhaps the simplest possibility. 
Gubser [27] has shown that coupling a neutral scalar field 
to the square of the Weyl tensor leads to a scalar hair 
on an asymptotically flat Schwarzschild black hole. This 
coupling is of interest here because the square of the Weyl 
tensor vanishes at the boundary (r — > oo) and hence does 
not affect the scaling dimension of the dual boundary 
operator. To this end, we modify the BGP Lagrangian 
(1) as follows 



J^ip A~y 



£■ 



BGP 



+ -rV^C^^p^C^'^P'^ i^.s^-^' - ^^) , (14) 



with C^i/pa the Weyl tensor. In the Appendix, we show 
that the above modification to the BGP Lagrangian, 
though not unique, offers a valid description of the neu- 
tral massive spin- two field (/3^jy, meaning that Lpp,^ still 
propagates causally in a fixed Einstein spacetime with the 
correct number of degrees of freedom. In what follows, we 
show that the equations of motion for ip^^j obtained from 
the new Lagrangian now admits non-trivial normalizable 
solutions (whose near horizon and asymptotic boundary 
behaviors are given by (12) and (13), respectively). 




FIG. 2. A plot of Bij verse A^ for f' = 1.0866 and m^ = 1/100. 
The values of Bij and Aij have been rescaled so that they are 



diniensionless. The gauge in (9) has been used to set ipx 



^xy- 



bound as before. Thus, regardless of the value of the cou- 
pling i"^ , the asymptotic AdS4 region of the background 
remains stable. 

As Figure 2 illustrates, there exists a normalizable so- 
lution with Aij ~ and Bij ^ 0. (We have used the 
gauge in (9) to set (pxx = '^xy) The holographic inter- 
pretation of this normalizable solution is that Oij has 
spontaneously condensed in the boundary theory. While 
it may be possible to obtain a normalizable solution for 
LPfj_i, using an alternative mechanism, we believe that the 
details of how the dual operator condenses in the bound- 
ary theory should be irrelevant to the Pomcranchuk in- 
stability as we discuss in the next section. 



D. Maxwell's Lagrangian 



C. Normalizable Solution 

From the Appendix, the relevant equation for deter- 
mining whether there exists a normalizable solution for 
^^iv is 



Since we will be considering the boundary theory at a 
finite chemical potential ^ for a U(l) charge, we need to 
introduce in the bulk an Abelian gauge field A^ with the 
Lagrangian 



C 



1 



M 



-gFnuF^^r 



(V^ - TO^ + fC-ySpaC^^P'^) ippy + 2RPp''^ipp„ = 0. (15) The equations of motion for A^, is then 



Note that the background metric must again satisfy the 
Einstein relation. Since we have taken the spacetime to 
be the Schwarzschild AdS4 black hole, the square of the 
Weyl tensor is Cp^p„C^"'P'' = 12(ro/r)6. One can easily 
show that, given the constraints presented in the Ap- 
pendix and the equation (15), it is again consistent to 
put (ppi, in the form given in (8). Notice that as shown 
above, the square of the Weyl tensor evaluated on the 
background vanishes as r — > oo. This then implies that 
the relationship between the asymptotic mass m of (pij 
and the UV dimension A of the dual operator Oij is still 
given by (11). Also, the mass of tpij in the asymptotic 
AdS4 region satisfies the same Breitenlohner-Freedman 



V^^^, = 0. 



(16) 



(17) 



We will only be interested in the case where At is non-zero 
and, moreover, depends only on the radial coordinate 
r. So, we take the ansatz Ap ~ At{r)5pQ. Given our 
assumption that Ap does not back react on the metric 
(7), the solution to (17) is easily obtained to be 



A = n{l-—\dt, 



(18) 



where the boundary conditions at the horizon and the 
asymptotic boundary have been fixed by demanding 
^(('^o) = and At{r -^ oo) = /x. 



III. FERMIONS 



In this section, we study the shape distortion of the 
Fermi surface as a result of the condensation of the 
boundary theory operator Oij. We start by introducing 
a bulk spinor field -0 with mass m^ and charge q which 
is dual to a boundary theory fermionic operator O^ with 
a UV scaling dimension A^ and charge q. Note that the 
bulk spinor f/^ is a four-component Dirac spinor while the 
dual operator O^ is a two-component Dirac spinor. In 
this paper, we take m^ € [0, |) and choose the so-called 
"conventional quantization" where the scaling dimension 
of the fermionic operator is related to the mass of the 



bulk spinor field through A^ — 



-m^. 



A. Fermion Lagrangian and Equation of Motion 

In our discussions below, the bulk spinor field "0 is 
treated as a probe where its backreaction on the back- 
ground metric, gauge field and neutral massive spin-two 
field is ignored. The Dirac Lagrangian, 



C 



Dirac 



= —^—giip{]p — m)'tp, 



(19) 



describes the bulk spinor field, where tp — i/i^ F- and 



]p = e^fF^ {df. 



1 , , ab 



^ab — iQ^fj.) a-iid e^ and i 



being 



respectively the vielbcin and the spin connection. Also, 
a,b, ■ ■ ■ = {t, X, y, r} denote the tangent space indices, 
F-, F-, F^, F- are the Dirac matrices satisfying the Clif- 
ford algebra {F",F^} = 2??''^ and F^f, = \[Ta,Tb]. The 
overall negative sign in (19) ensures that the bulk ac- 
tion, once holographically renormalized, gives rise to a 
positive spectral density for the fermionic operator O^. 
In (19), g^ij and yl^ are the background metric and the 
background gauge field whose expressions are given in (7) 
and (18), respectively. 

In the probe limit that we are considering here the 
Dirac Lagrangian (19) is not capable of holographically 
realizing how the spectral function of the boundary the- 
ory fermionic operator O^ is affected in the presence the 
order parameter {Oij). To do so, one has to explicitly 
introduce bulk couplings between the spinor field ip and 
the neutral massive spin- two field (p^^. A similar problem 
arises in the context of holographic d-wave superconduc- 
tivity [22, 23] although in that case one has to deal with 
coupling the bulk spinor to a charged massive spin-two 
field. Since our goal here is to capture the leading order 
effects on the retarded two-point function of the fermionic 
operator due to the presence of the order parameter, we 
consider those couplings in the bulk which are quadratic 
in the spinor field and have relatively low mass dimen- 
sion. Furthermore, among the aforementioned couplings, 
we will only be interested in those which will potentially 
give rise to anisotropic features in the fermion spectral 
function. In other words, we ignore the terms which 
would just modify the mass of the spinor field. Leav- 



ing aside those terms which would vanish once evaluated 
on the background, we focus on the leading-order"^ non- 
trivial coupling (which can potentially contribute asym- 
metric features to the fermionic spectral function) 



>Cint = -i>^V--9fiiA^^D'''il), 



(20) 



with A real. Thus, the bulk fermion Lagrangian that we 

study is C^ == ^Coirac + >Cint- 

The Dirac equation following from C^ takes the form 



(0 - m + \^^,T^D'')ij = 0. 



(21) 



To solve the above equation, we go to momentum space 
by Fourier transforming tp{r,x^) ^ e'^^'^^{r,k^)^ where 
fc^ — {u!,k). It is expedient to choose a basis for the 
Dirac matrices as follows 



pa^ 



-0-3 
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f2 
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(22) 



It is also convenient to remove the spin connection 
from the Dirac operator Ip by rescaling ip according 
to -0 — > (— g5^'')^^/^-0. Splitting the new (rescaled) 
spinor ip according to ip'^ — (ipi,ip2), where tpi{r;uj,k) 
and il)2{r;Lo,k) are two-component spinors, and given 
the ansatz for the background form of ip^^, where the 
only non- vanishing components are ipxx — ~'^yy £^nd 
'fxy ~ 'i^yx, the equation (21) results in the following 
two coupled differential equations for ipi and -02 

(A- - i(T2y/9^k) tpi + i(J2^/9^~ky^2 = 0, (23) 
A- + i(y2^/g^K) il>2 + i<J2Vg^kyi;i = 0, (24) 



where we have defined 



Dr = -(Tsy/^dr + CTl V^ (^ + 1^2) 

kx^kx{l + Xg^^'^iPxx) + >^g''''Vxyky, 

ky^ky{l- Xg'^'^Lpxx) + A.g^^(^a:j,fca;- 



TO^, 



(25) 
(26) 
(27) 



For A = 0, or when ip^x = Vxy = Oj the Dirac equa- 
tion (21) has been studied extensively in the literature 
(although mainly on the Reissncr-Nordstrom AdS4 black 
hole background), following the work of [28-31]. These 
studies show that in these cases there are symmetrical 
Fermi surfaces in the boundary theory whose underlying 
excitations can either be Fermi or non-Fermi liquid-like. 
On the other hand, when ip^x develops a normalizable 



^ There is also a coupling of the form As-y^— g ip^vi> T^V^ D^ ip that 
breaks parity in the boundary theory when the operator dual to 
ip^^ condenses. Wc will not consider such a coupling in this 
paper. 



mode in the bulk, k^ and ky will differ thereby giving 
rise to asymmetrical features in the fermionic correlators 
of the boundary theory. This can be shown explicitly. 

To compute the retarded correlator of the boundary 
theory operator O^, one needs to solve the equations (23) 
and (24) with in-falling boundary condition at the hori- 
zon [32] and read off the source and the expectation value 
of O^ from the asymptotic expansion of "00. (a = 1, 2) 
following the prescription of [33]. Indeed, choosing in- 
falling boundary conditions for f/^Q, near the horizon, the 
leading-order asymptotic [r — ^ oo) behavior of i}}^ takes 
the form ipair -^ oo) — {Bar~^^'t' .Ao.r"^''')'^ , where the 
two-component spinor A = (Ai, ^2)"^ sources the bound- 
ary theory fermionic operator 0^, while B = {Bi,B2)^ 
gives the vev of the operator. The spinor B is related to 
the spinor A through B = SA from which the retarded 
Green function of the operator O^ take the form [33] 



Gfl(c^,fc) = -z57* 



(28) 



Note that in our basis of gamma matrices 7* = zcti . Up 
to a numerical constant, the spectral function of the op- 
erator O^p is given by p(a;, k) — TrImGij(aj, k). 



B. Fermi Surfaces and Low-Energy Excitations 

To obtain the retarded Green function (28) of the 
fermionic operator 0,0, one should numerically solve the 
equations (23) and (24), or the flow equations obtained 
from them, with in-falling boundary condition at the 
horizon. We do this computation at finite temperature 
mainly because our background, in which there is a nor- 
malizable solution for the neutral massive spin-two field, 
is not, strictly speaking, valid at zero temperature. The 
Fenii momentum fcp appears as a pole in Re Gj^{uj, fc), al- 
though at finite temperature such poles, instead of being 
the delta function, become broadened. 

Carrying out the numerical computation, we find mul- 
tiple Fermi surfaces for the parameters chosen. To 
demonstrate proof of concept, it is sufficient to focus on 
the first Fermi surface, where \k-p\ (measured compared 
to the effective chemical potential g/i) has the smallest 
value. In generating the plots shown in this section, we 
first used the gauge (9) to set ip^x = fxy and then nu- 
merically solved the Dirac equations for a very small but 
non-zero frequency such that the delta functions at the 
Fermi surface are broadened. 

First, consider the case where the background value of 
ftiuiT) identically vanishes in the bulk. In other words 
the boundary theory is in the symmetry-unbroken phase 
(Oij) — 0. Not surprisingly, the Fermi surface obtained 
in this phase is isotropic, as is evident from Figure 3(a). 
To investigate the nature of the excitations, we com- 
puted the quasiparticle dispersion relation by focusing 
on the peak in the spectral function. In a Fermi liquid, 
uj{k) = k^ with k± = k — kp and z — 1. Figure 4 (see 
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FIG. 3. Density plots of the spectral function of the fermionic 
operator O^ in the boundary field theory, a) Spectral density when 
the boundary field theory is in the unbroken phase b) Spectral 
density in the broken phase showing an elliptic Fermi surface where 
we set A = —0.4, c) spectral density in the broken phase where 
A = —0.95 showing a suppression of the spectral weight and d) 
spectral density, also in the broken phase, when A = — 1.4. 



open circles) reveals that in the unbroken phase of the 
boundary theory where the Fermi surface is circular, the 
dynamical exponent z — 0.64 which indicates that the 
low-lying excitations form a non-Fermi liquid. Impor- 
tantly then, any subsequent breaking of rotational sym- 
metry of the Fermi surface will be from a non-Fermi liq- 
uid state, one of the key hurdles in describing the strong 
correlations in the experimental systems. 

Focusing now on the interesting phase where the neu- 
tral symmetric traceless operator Oij has spontaneously 
condensed, we find three significant features in the spec- 
tral function of the fermionic operator O^. First, the 
rotational symmetry of the Fermi surface is broken. As 
we alluded to earlier, such behavior is expected given the 
form of kx and ky at each r-slicing in the bulk. Second, 
the spectral weight is angularly dependent. For interme- 
diate values of the coupling A, a gap-like feature opens 
at the two end points of the major axis. This highly 
non-trivial behavior is not a generic feature of traditional 
treatments of the Pomeranchuk instability [8] but can be 
understood within a WKB approximation to the spec- 
tral function as outlined in the next subsection. The 
diminished spectral density reappears once A exceeds a 
critical value as shown in Figure 3(d). Finally, the excita- 
tions near the Fermi surface remain non-Fermi liquid-like. 
The solid circles and triangles in Figure 4 show that the 
dispersions along kx and ky deviate from linearity with 
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FIG. 4. The dispersion relation, Lo{k) oc fc^, of quasi-particles 
around fcp plotted on a log-log scale. The three cases correspond 
to a) open circles for the unbroken phase of the boundary theory 
with a circular Fermi surface where we find z = 0.64, b) solid 
circles for the dispersion relation (in the broken phase) along x- 
axis where we find z = 0.74 with A = —0.4 and c) solid triangles 
for the dispersion relation along y-axis (in the broken phase) where 
2 = 0.84 with A = -0.4. 



exponents oi z — 0.74 and z — 0.84, respectively. These 
values are, of course, parameter dependent. 



C. WKB Analysis of Spectral Function 

In this section we employ a WKB approximation to 
uncloak uncloak the origin of the momentum-dependent 
spectral weight in Figure 3. Such an approximation has 
been used recently in [34] for the semi-analytic analysis 
of spectral functions in the electron star background [35] . 
Setting kx — k cos 9 and ky = k sin 0, we define new 
momentum coordinates 



K^ 



1.2 



k'y^O, 



(29) 






which are written in the so-called "lab" frame. Note that 
the expressions for kx and ky have been defined in (26) 
and (27), respectively. With these new momenta and for 
a particular 9 and fc, the original Dirac equations (23) 
and (24) give rise to two decoupled equations 



(30) 
(31) 



Any anisotropy that arises now in the spectral function 
must be tied to k'^- To confirm this, we construct a con- 
tour plot of k'x at the horizon where the anisotropy is 
largest. The plots in Figure 5 show clearly that for a 
fixed value of k, k'^ is minimum at 9 = 7r/8 and Ott/S 
which is consistent with the values of 9 where the spec- 
tral density is a minimum. 

Indeed, one can develop a better understanding on how 
k'x affects the spectral function by performing a WKB 
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FIG. 5. Contour plot of fc^ for a) A = 0, b) A = -0.4, c) A = -0.95 
and d) A = —1.4. The plots have been generated with m = 0, oj = 
and g/i = 17. 



analysis of the Dirac equation. In order to make the 
notation less cluttered, we will replace fc^ by k in what 
follows. In order to apply the WKB approximation, we 
will take the limit where fc, to, uj and q are large, but their 
ratio remains constant. To this end, it is convenient to 
rescale all of these quantities with a factor 7 such that 
k = jk, m^ — jrh, q = jq and w = 7W, where 7 is taken 
to be large (compared to one). 

Central to the WKB approximation is the construction 
of a Schrodinger-like equation for each of the components 
of the two-component spinor 



V'l 



$1 

$2 



The same treatment applies equally to 'ip2 with the re- 
placement of fc — > — /c. To leading order in 7, the second- 
order differential equation for <I>i takes on the form 



1 = 7^5"- |m^ + g" [io + m (1 - rQ/r)f 



$'/ = 



with 



$„ 




-\/^(w + g^*) 



(32) 



(33) 



In the context of a Schrodinger equation, the right-hand 
side of equation (32) can be regarded as the zero-energy 
potential 

V = j^grr{m^ + g" [^ + A*g(l - ro/rf + 5"^}. (34) 




horizon 



FIG. 6. Potential V in the Schrodinger equation as a function of 
the radial coordinate r. The classical turning points as measured 
from the horizon are denoted by ra, r2, and ri. 



It is the turning points of this potential that governs the 
physics of the WKB approximation. In the near-horizon 
region, the potential is approximated by 



V{r 



ro) = 



3rj](r-ro) 



L^P + m^rl - uj^ 



3(r-ro) 



(35) 



If to is non-zero, the near horizon potential is dominated 
by the term proportional to w^ and is negative for real 
positive Lo. Close to the boundary, the leading order term 
of the potential 



V{r 



2 7-2^2 



j^L^m 



(36) 



is always positive. Between these two regions, there can 
be at most three turning points which we denote by ri, 
r2 and r^, with ri being the turning point closest to the 
boundary as depicted in Figure 6. For the parameter 
range of interest here, all three turning points will enter 
the matching conditions as we describe below. 

What remains is a simple matching procedure to de- 
termine the form of $i in each of the regions that bracket 
the turning points. This procedure, detailed previously 
in [34], utilizes the functions 



rr2 

X ^ drW + \og2, 

J rs 




(37) 
(38) 



which directly enter the retarded Green function 



GR{oj,k) ex — lim r 

Z r— foe 



2mL, 



where 



cxp(-2 / dr'W] (39) 



(40) 



cosh(X + iY) + sinh(X - iY) 
cosh(X - iY) - ismh{X + iY) ' 



The poles of the retarded Green function only come from 
Q and are given by y = irn — ie~'^-^ -!-••• where we have 
included the imaginary part up to the leading order only. 
Since we are interested in the non-analytic dependence of 
the Green function on frequency and momentum, we will 
focus on the behavior in the vicinity of the pole where Q 
can be written as 



Q = 



^-^ Y{uj, k)-Trn + ie-2X((^,fc) ' 



(41) 



p , we can expand Y{u;,k) as 



Close to uj = and k ~ k^' 

As a result of this expansion, the retarded Green function 
takes the form 



GR{uj,k) ex ^ 



u(")^ 



E 






(42) 



where 



[w + v,,{k-kf^)Y+cle-^x 

" 9.r(o,4"V 



a.r(o,4")) 



(43) 



j 

dr^V{Q,kP), 



which is essentially the equation (48) of [34]. 

Close to the Fermi surface, where the first term in 
the denominator of the second equation in (42)) can be 
dropped^, the spectral function j4(w, k) can be extracted 
from the imaginary part of the retarded Green function. 
One then obtains 



A{lo^ k) ex e 



-2a„+2X 



(44) 



We can now understand the anisotropy of spectral func- 
tion from the dependence of a„ and X on the area en- 
closed by the potential integrated between the turning 
points r2 and r^^ and between ri and oo which will de- 
termine a„. In fact, for ttz = 0, ri — > oo, a„ therefore 
vanishes and the angular dependence of the spectral func- 
tion will be dominated by X alone. Close to the horizon. 



Note that in equation (34) we have written the potential in terms 
of the hatted quantities Co, k, however, in doing the expansion of 
Y{io, k), the factors of 7 cancel and we can drop the hat on ui 
and k. 
We have assumed that in the limit uj —> 0, io vanishes faster than 



-2X(l. 



'*' which is true in the range of parameters in which we 



are interested. 

This will define the quantity X. 
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FIG. 7. Effective WKB potential for A = -0.4, m = at the first 
Fermi surface. The solid line corresponds to 6 = and the dashed 
line to e = 7r/8. 



the first turning point of tlic potential, r3, is rouglily 
given by 



rs « ro + 



w^L^ro 



3[m2r2 + fc(ro)2i2] 



(45) 



From tliis we can see tlrat in the limit w — >■ 0, that is, 
close to the Fermi surface, r^, — > tq. 

In our analysis of the spectral function, where we 
have set ?7i = 0, only X determines the spectral func- 
tion. To illustrate the angular dependence, we plot the 
near-horizon effective WKB Schrodinger potential for 
A = —0.4, TO = at first Fermi surface (n = 0) in Figure 
7 (corresponding to Figure 3(b)). The dashed line cor- 
responds to the potential along 9 = 7r/8, while the solid 
line corresponds to = 0. It is now obvious that r2 is 
larger for 6* = than for 6 — tt/S and the potential is 
always larger for ^ = than 9 — tt/S in the range of 
r shown. Hence X and spectral weight are larger along 
9 = than 9 = 7r/8, thereby explaining our key finding 
that the spectral weight is angularly dependent. In prin- 
ciple, the WKB approximation works only for large n, 
nevertheless, we see that the result qualitatively agrees 
with our numerical results even for small n. 



IV. DISCUSSION 

We have shown here how holography can be used to 
model a shape distortion of the underlying Fermi surface 
in a non- Fermi liquid through the condensation of a neu- 
tral symmetric traceless operator which is dual to a neu- 
tral massive spin-two field in the bulk. An open challeng- 
ing problem remains: Is the condensation of the bound- 
ary theory operator Oij possible starting from a bulk 
geometry which is not an Einstein manifold? In other 
words, what is the backreacted background that allows 
for the ghost-free and causal propagation of a neutral 
massive spin- two field. These extensions are extremely 



desirable since they would not only enable an analysis at 
zero-temperature but would also cure at finite tempera- 
ture some of the peculiar thermodynamical properties of 
the conduced phase in the probe limit. In addition, is 
the angular dependent spectral weight a generic feature 
of holographic non-Fermi liquids with (partially) broken 
rotational symmetry? If so, then holography would have 
provided a key ingredient missing from most condensed 
matter analyses of the Pomeranchuk instability, namely 
a vanishing spectral weight for a finite range of momenta 
and hence a pseudogap as a function of frequency. Also, 
how does the resistivity tensor behave in the two different 
directions? In our probe analysis, where the backreaction 
of the massive spin-two field on the metric is ignored, 
the resistivity tensor is not sensitive to the anisotropy 
of the system when the boundary theory is in the bro- 
ken phase. Nevertheless, one might be able to see the 
anisotropic contribution of the fcrmions to the resistivity 
through a fermion loop computation in the bulk similar 
to the analysis of [36]. 
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APPENDIX 

In this appendix, we demonstrate that adding a cou- 
pling of the form fC^spaC''^'""i^^^u^'"'-'p'^) to the BGP 
Lagrangian, given that the background spacetime is an 
Einstein manifold, does not induce any extra degrees of 
freedom for the neutral massive spin-two field tp^i, nor vi- 
olate its causal propagation. Our discussion below closely 
follows that of [19], where it has been shown that neu- 
tral massive spin- two field in an Einstein spacetime, such 
as the AdSd+i Schwarzschild black hole, has the correct 
number of degrees of freedom and propagates causally. In 
our discussion below, we keep the number of dimensions 
of the bulk spacetime (denoted hy d+1) arbitrary. 

Consider the following Lagrangian [20] 



C = 



+ 2a^Rip^,,ip^''' + 2a2Rp^ + 2a^R''^''P^^,^xp 
+ 2a4i?^VMA'/'^ + Sasi?^ Vp.^^} , (46) 



10 



where ip = ip^t. The above Lagrangian is the most gen- 
eral two-derivative action for a neutral massive spin-two 
field ipfj^i, (up to the quadratic order in f>^^, though) in a 
curved background. The coefficients in the first two rows 
of (46) are fixed because in the flat spacetimc limit, one 
demands the Lagrangian to go over to the Fierz-Pauli La- 
grangian [18], which is known to have the correct number 
of causally-propagating degrees of freedom for a neutral 
massive spin-two field. Indeed, in order for a massive 
spin- two field (represented by a symmetric tensor (p^v) 
to have the correct number of propagating degrees of 
freedom in a (d -I- l)-diincnsional spacctime, one needs 
2{d + 2) constraint equations to kill the extra degrees of 
freedom. 

Now, add to the above Lagrangian an additional cou- 
pling of the form £'^C'^{ipf^i,ip^'^ — ip'^) and define the La- 
grangian 



C^=C+ l^C^ [p^,^^' 



(47) 



the covariant derivative of the equations of motion 



(50) 



Again, neither V^E^y nor V^ \C'^ {ip f^^ip^" — </?^)] in- 
volves two time-derivatives of the spin-two field. Thus, 
the equations (50) give us d -I- 1 secondary constraints. 
Following [19], the equations (50) should contain (p^Q 
through a rank d matrix in order for their derivatives 
to define d accelerations (pi^: 



^MO,^ 



E„ 



Ap'm +B^p'jo 4 
^ ^Ciip'oQ +T>fip'jo - 
A B^ 



(51) 



rank $"" 



rank 



C,. Pf 



The explicit form for $5^ when £ ~ appears in [19]. 
additional non- vanishing contributions to A, B-' , Ci 
2?* are given by 



SB^ = -tC^g- 



2„j0 



5V) ^ eC'^g^ 



^5). 



The 
and 

(52) 



where, in order to make the expressions less cluttered, 
we defined C^ = C^gp^C^^P" . For ^ = 0, the analy- 
sis has been performed in [19], where, in order to get 
rid of the extra degrees of freedom of Pp,ij, one obtains 
a one-parameter family of solutions for the coefficients 
Oi,--- ,05 (labeled by an arbitrary real parameter ^) 
which can be represented as 



ai 



1' 



a2 



1-2^ 

2d + 2' 



as = 04 = as = 0, (48) 



provided that the background metric g^^ satisfies the 



Einstein relation Rr^, 



2Kg^y/{d 



2(^ - l)/{d + 1). Defining m^ 



2(1 



and rriQ 7^ 

mi{d+i) 



and given the above solution for a^'s, the Lagrangian (46) 
takes the form given in the expression (1), the so-called 
BGP Lagrangian. 

For £ 7^ 0, one can easily obtain the constraint equa- 
tions following [19]. Basically, constraints come from the 
equations of motion and their derivatives which do not 
contain two time derivatives of spin-two field p^u- The 
equations of motion read 



E„ 



E^,^ + tC^ {pf,^ - gpyip) = 0, 



(49) 







where i?^^ denotes the equations of motion for the 
case. The explicit form of _E^^ is given in [19]. Since 
EpQ does not involve second time-derivatives of p^i, , and 
obviously the new term in (49) does not involve any sec- 
ond time-derivative of (p^i, either, the equations E^q ~ 
actually give us d -(- 1 (primary) constraints. 

In order to obtain the secondary constraints, we take 



Hence, in order for the rank of $5^ to be d, we require 
ml - eC^ + 2(1 - C)Rl{d + 1) 7^ 0, where ^ is arbitary 
parameter. We now have 2d -I- 2 constraints. 

One of the remaining constraints can be obtained from 
a linear combination of the equations of motion and the 
primary and secondary constraints: 



ml 



PC^ ~ 



d-1 
d-1 



E^^ 



yp^VE^ 



(d + l)(d-l) ^ 



1 



2(1-0 
d+1 

-2Cd 



R 



tc"- 



d + 1 

2^2(v''C2)(V^</.^, 



R+{ml-tC'-)d 



(53) 



^M 



+ ^2(vMv-c2)(^^,-g^,(^)«0. 

This equation can be used as a constraint since it does 
not contain py^^,. It contains pp_i, which can be removed 
by the equations S/^E^^i, = 0. After such a removal, the 
time-derivative of (53) does not contain (/3^^ and can be 
used as the last constraint. We then have in total 2d -I- 4 
constraints and hence the correct number of propagating 
degrees of freedom. 

Following the discussion in [19], it can be seen easily 
that the characteristic matrix remains unchanged for £ 7^ 
0, due to the fact that the Weyl squared coupling does 
not affect any of the terms with two time-derivatives. So, 
the equations remain hyperbolic and causal even for the 
£ ^ Q case. The resulting equations of motion are 

(V^ -m^+ ec'') ^p, + 2Rp;yP^p, = 0. (54) 

where, as we have defined before, m'^ = ml + 2(1 — 



e)i?/(d+i). 
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